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TECHNICAL NOTE 3165 

USE OF A HOT-WIRE ANEMCMETER IN SHOCK-TUBE INVESTIGATIONS 
By Daxshan Singh Dosanjh 


SUMMARY 


Thp use of the hot-wire aaemometer as an instrument for experimen- 
tal investigations of the transient flow phenomena in a shock tube is 
examined. The response of a hot-wire to a transient step-function type 
of change in flow conditions (e.g., shock waves) is experimentally stud- 
ied and a criterion for predicting the right compensation is postulated 
and experimentally verified. The advantages and limitations of this 
somewhat special and, xmtil now, unattempted use of the hot-wire anemom- 
eter are stressed. 

Besides using the hot-wire anemometer for the actual wave phenomena 
and for flow studies, it has been very successfully used as a timing 
and/or triggering device in shock-tube work. It is demonstrated both 
from theoretical considerations as well as experimental observations that 
the hot-wire anemometer is satisfactory for ti m ing very weak traveling 
shock waves, while it can be used for similar purposes for comparatively 
stronger shocks. 


INTRODUCTION 


The shock tube as an instrument of research has been gainfully 
employed in the studies of one -dimensional unsteady flow of gases and 
the transient wave phenomena by a large number of workers both in this 
country and abroad. The shock tube has served as an inexpensive and 
simple sovirce of transient flows of quite heterogeneous thermodynamical 
properties. The short-diiration flows with pressures, temperatures, and 
densities extraordinarily different from the normal atmospheric condi- 
tions are easily made accessible. For instance, one can easily produce, 
for a short duration, temperat\ires ccanparable with surface solar 
temperatxares . 

Optical techniques^ (spark shadowgraphs, schlieren, interferometer, 
spectroscopes, etc.) have been the main tools employed in shock- tube 

^References 1 and 2 ; also private conversation with Research Group 
at Cornell University, Department of Aeronautics, about some spectro- 
scopic work being done in the region behind traveling shock waves. 
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research. These techniques have been generally very useful in studying 
the traveling shocks, their reflection and diffraction, and other allied 
phenomena. But if there were some fluctuations in the main flow behind 
these thermodynamic discontinuities (shock waves, contact surfaces, 
rarefaction waves, etc.), then these instruments shed only a limited 
amount of quantitative light on the phenomena (ref. 5)* For instance, 
it is well known that because of the breaking characteristics of cello- 
phane (the usual diaphragm material) the contact surface flow is rather 
turbulent. Qualitative confirmation of the presence of turbulence in 
this flow is the best that optical techniques have been able to deliver 
so far. It was to explore the possible eventual use of the hot-wire 
anemometer for the study of transient mean and/or fluctuating flows and 
the wave phenomena involved in shock tubes that this present work was 
undertaken. 

The hot-wire is a standard instrument for studying the nature of 
subsonic and supersonic turbulent flow fields. Almost all of the devel- 
opment and worth-while analysis of this commendable technique are based 
on the assumption of small fluctuations as compared with the mean flow 
in which the hot-wire is ox>erated. Ordinarily shock tubes are operated 
such that there is no mean flow in the system to start with. Then all 
of a sudden large flow fields are established for a very short duration. 
If the hot-wire is mounted in the middle of the cross section of the tube 
at some convenient location and the shock tube is operated, the response 
of the hot-wire will give the time history of the flow passing over it. 

So before one is able to study the fluctuations (artificially introduced 
or naturally present) in the flows behind the principal waves generated 
in the shock tube, it is necessary to gather some basic information about 
the response of this instrument to these transient wave fronts and the 
flow regions behind them. 

The hot-wire anemometer and the associated electronic equipment are 
usually tested on the basis of the validity of certain hot-wire response 
equations (refs, k, 5, and 6). The transient response of the hot-wire, 
when proper electronic compensation is applied, is usually assumed to be 
identical for transients in flow conditions or in heating current. 

(Ref. 7 Betchov seems to show that, because of thermal propagation, 
the two may be somewhat different.) This plausible assumption was first 
experimentally verified by the author and others (ref. 8). Since then 
the technique of experimentation has been considerably improved and most 
of the resialts arrived at in the original report have been rechecked and 
extended. Main conclusions included in the present report are based on 
data obtained mostly with the experimental arrangement as described here. 

Known fluctuations in velocity and temperature are difficult to pro- 
duce. Wires have been shaken sinusoidally up to I 50 to 200 cps in moving 
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air streams, but the frequency was limited because of the mechanical 
difficulties.^ It is believed that the traveling shock wave furnishes 
an ideal flow field for stiidyiiig the hot-wire response to a time -step- 
function change in flow conditions. When a traveling shock wave passes 
over the hot-wire, the latter experiences a discontinuous step-function 
variation of the mass flow, pressvire, density, stagnation, and ambient 
temperature, and the flow is uniform and constant thereafter for a few 
milliseconds until the contact surface or some other wave front (reflected 
shock or rarefaction wave, etc.) reaches the hot-wire. 


When a step-function type of transient (either in cvirrent or flow 
condition) is applied to the hot-wire, because of the thermal lag the 
response of the hot-wire is exponential (ref. 6). This can be electron- 
ically compensated such that the exponential is restored to almost a step 
function (ref. 8 and fig. 1). The hot-wire anemometer thus can be used 
to study the transient wave fronts and fluctuations in the flow fields 
in a shock tube if one knows how to predict the right compensation values 
(for constant-current operation) . The criterion for predicting the com- 
pensations has been sviggested and verified within an acceptable experi- 
mental agreement. 


Ultimate time resolution of the hot-wire and the associated equip- 
ment is usually estimated electrically (square-wave technique, ref. 4). 
Hot-wire response to a shock wave offers an alternate way of experimen- 
tally determining this \iltimate time resolution by recording its response 
to a controlled step-function type of change in flow conditions. This 
has been done for the equljanent used in this work (see also ref. 8). 

In shock-tube studies, one of the most Important Instrumentation 
problems is the controlled timing of the trainsient phenomena. Usually 
the primary shock wave is used as the basis for the design and operation 
of a transducer which changes some property of the shock front into a 
suitable electric pulse. This pulse, after sharpening and amplification, 
is utilized for triggering and/or timing purposes. Since the exponential 
response of the hot-wire to any sudden changes in flow conditions can be 
sharpened either by ccanpensating amplifiers or by differentiating the 
signal, the hot-wire anemometer has been successfully used to time various 
transient phencanena in a shock tube. 


Most of the experimental work in shock tubes does not include any 
systematic studies at the lower end of the spectnam of shock strength. 
This, by and large, is due to the limitations of the existing timing 
devices. When the hot-wire is used at elevated temperatures, it is pre- 
dominantly mass flow sensitive. This characteristic of hot-wires has 


p 

^Unpublished i>ersonal communication regarding work done on low fre- 
quency shaking of hot-wires at the National Bureau of Standards in the 
early 1950 's. 
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been successfully used to time very weak shocks (drift velocity of a few 
feet per second) nearing virtually a sound speed. The experimental 
triggering anti timing system actiially used is described in the text. 

(For details see ref. 9 *) 

It is experimentally demonstrated over a limited range of shock 
strengths that the hot-wire anemometer can be used to measure transient 
temperature jumps. Also the mass-flow and temperature jumps which occur 
simultaneously in case of a shock wave can be separated. This is achieved 
by using two different heating currents thro\igh the same wire, with the 
same shock strength used twice consecutively (ref, 5 )- The difference 
between the observed results and the expected ones is discussed. Hot- 
wires have also been used in detecting and timing rarefaction waves. 

Since this technique has further potentialities, the discussion 
included about the response of the hot-wire in various regions of the 
shock-tube flows at times exceeds the reinge and regions actually inves- 
tigated experimentally. This gro\md work may facilitate the work of those 
who may choose to use this technique to extend the scope of the present 
investigation . 

The author is indebted to Drs. Francis H. Clauser and Leslie S. 

G. Kovd,sznay for their encouragement and many helpful suggestions. 

Mr. Richard Swartley's help with the experimental work is gratefully 
acknowledged . 

This work was supported imder the combined sponsorship of the Bureau 
of Ordnance, Department of the Navy, and the Department of the Air Force. 
Since April 1, 1955^ "the sponsorship has been entirely by the Air Force. 
The author appreciates this much needed financial assistance. 

The present report is based on Part I of a dissertation (ref. 10) 
submitted by the author to the Faculty of Philosophy of The Johns Hopkins 
University in conformity with the requirements for the degree of Doctor 
of Philosophy. It has been made available to the National Advisory 
Committee for Aeronautics for publication because of its general interest. 


SYMBOLS 

Hot-Wire Symbols 


av(t) ^ - cte(TY - Tg) 

«e 

a,(t) = . Oa(T„ - Te) 
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C 

d 

E(Tw) 

H 

I 


k 

I 


M* 

n = -^ 

aRjp 

Re(t) 

Rf 


heat capacity of wire, ergs/deg 
diameter of wire 
stored energy in wire, CT^ 
voltage fluctuation or jump 
heat loss from wire per unit time 
heating current through hot-wire 

characteristic current of wire (a calibration constant used in 
range where King's equation of heat loss is assumed valid) 

heat conductivity of gas 

length of wire 

time constant of wire for small fluctuations (linearized theory) 

equilibrium resistance at T^(t) 
reference resistance at 

wire resistance at room temperature without any current 
through it 


Ry(t) wire resistance at Ty(t) 

R^(t) wire resistance at T^(t) 

Tg(t) equilibrium temperat\ire attained by wire if unheated 

Tf reference temperature (usually at 275^ K) 

Ty(t) virtual wire temperature that would be attained in absence 
of thermal lag 

T^(t) instantaneous temperature of hot-wire (average along wire) 

t time 

Uq characteristic velocity of wire (a calibration constant of 

wire used in rajoge where King's equation of heat loss is 
assumed valid) 
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W heat generated in wire per \mit time 

a temperature coefficient of resistivity at reference temperature 

a 

= 

1 + a(Tg - Tf ) 

H viscosity of gas 

a 


Shock-Tube Parameters 

local speed of sound 

^vc ~ ^oc 
^oc 


„ Rv2 - ^e2 _ Rv2 - 1^02 ,4. i. _ 

av2^t) = = xf Rg2 - 


Cy 

K 

M 

M, 


2 ,c 

Mou*(t) 

Nu 

P 

Pr 

Pr 

R 

Re 


-p 'o 2 

^e 2 ^o 2 

specific heat at constant pressure 
specific heat at constant volume 
Knudsen momher, A/d 

Mach number of flow corresponding to region indicated by 
subscripts, u/a 

*(t) time constant for large step function in flow conditions when 


wire is operated "hot” 

time constnat for large step function when wire is operated 
"cold" at velocity u 

Nusselt number, H/:jtlk(Tv - Te) 

pressure ratio across diaphragm, Pj^/P]_ 

Prandtl number, M-Cp^k 

pressure in compression chamber 

resistance of hot-wire at corresponding T 

Reynolds number , pud / [x 
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Ry(t) corresponding resistance at T^(t) 

S measure of shock strength, P 2 /P 1 > 1 

T,p^p temperature, density, and pressure, respectively 


Tv(t) 


= Constant for t < 0 

operating equilihriijm temperature of 
T^2 = Constant for t = 0 wire after wave front under study 

has passed over wire 


U 


s 


Us I 

u 



primary shock speed 
reflected shock speed 
velocity of flow 
coordinate along axis of tube 
l.k for air 


A molecular mean free path 

pu mass flow behind wave under study 

7 + 1 ^ « 

a = = 6 for air 

7-1 

Subscripts: 

c contact -surface region ending at foot of rarefaction wave and 

corresponding parameters 

e equilibrium value of parameter 

o stagnation parameters 

R region ahead of head of centered rarefaction wave; before 

bursting of diaphragm, compression chamber as a whole 

r region between head and tail of centered rarefaction wave 

V "virtual” parameters, that is, value of various parameters if 

there had been no lag due to thermal capacity of wire or if 
fluctuations were infinitely slow 
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region and corresponding parameters in front of primary shock 
wave; before shock wave is generated, conditions in expansion 
chamber as a whole 


2 


region from primary shock wave to start of contact-surface 
region and corresponding parameters 


5 


region between reflected shock and reflecting wall and corre- 
sponding parameters 


BASIC WAVE MODEL IN SHOCK TUBE 


A shock tube usxially consists of a long metallic pipe or box sepa- 
rated into two sections by a membrane across which a pressure differen- 
tial can be established. The high-press\ire side is called the ccmpres- 
sion chamber and the low-pressure side, the expansion chamber. When a 
desired pressure ratio is attained, the membrane is ruptured, generating 
a pressure wave propagating down the expansion chamber. Within a short 
distance from the membrane (a few diameters of the tube), the pressiire 
wave steepens into a sharp pressure discontinuity (shock wave). The 
simple wave model and the associated regions of flow assumed in the sub- 
sequent analysis of shock-tube phencmena and the behavior of hot-wires 
are reproduced in figure 2 . Mainly it consists of: 

(a) A centered rarefaction wave moving to the left in the compres- 
sion chamber. 

(b) A plane shock wave (that will be referred to as the primary 
shock wave) moving at constant speed to the right in the expansion 
chamber. 

(c) A contact surface (the boundary between masses of the gases 
previously on the two sides of the membrane) following the shock wave 
into the expansion chamber. 

(d) A reflected shock from the closed far end of the expansion 
chamber, traveling upstream. If the end is open, the reflected wave 
will be a rarefaction wave. 

Some of the other wave fronts in the shock tube are overlooked here 
because they do not play a prominent role in the scope of this work. 

The usual simplifying assumptions in the analysis of the various 
waves arid the associated flow phenomena are as follows: 

(1) Perfect gas with 7 constant. 
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(2) The membrane, or what is commonly called a diaphragm, is a plane 
boundary that is instantaneously removable without loss of energy. 

(3) The flow in various regions of the shock tube is isentropic and 
one dimensional. 

(4) There is no diffusion across the contact siirface. 

(5) The regions labeled R, 1, 2, and c as shown in figure 2 are 
assumed to be constant. 

The operation of both the compression and the expansion chamber with 
air in them will be Indicated by the notation Alr/Alr. 

The quantitative values of the various thermodynamic parameters 
(mass flow, temperature, etc.) in the various regions of interest are 
required for theif subsequent application in the determination of the 
hot-wire response. However, only the relevant relations without any 
extensive derivations will be included. (For details see refs. 1 and 10.) 


State 2 is the region behind the primary shock wave and the contact 
surface (fig. 2). All the following relations refer to the coordinate 
system in which the medium ahead of the shock front is at rest. 

Let the shock strength 


Flow Region Behind Primary Shock Wave (State 2) 


^2 _ Pressure behind shock wave 


Pressiire in front of shock wave 


and the primary shock speed 



( 1 ) 


where 


o 


7 + 1 
7 - 1 


A plot of Ug versus S is given in fig\ire 3* 
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The mass flow behind the shock wave in region 2 is given as 

P2^2 ^ (a - 1)(S - 1) l l + ~aS (2) 

Pl^l (T + S ^ a + 1 

A plot of Pp^p/Pi^i versus S is given in figure 4. 

The temperatirre ratio is 

^ = S(g + S) (5) 

Ti 1 + oS 

Since later^ the hOt-wire will be assumed to respond to stagnation temper- 
atures, it is desirable to obtain the corresponding relations (ref. 8 ). 
For the flow behind the primary shock wave, the energy equation is 


CpT2 + i U2^ = CpTo2 


or 


E °2 = 1 + 2 IJ 1 I 

T2 2 [^2/ 


One can get 

^02 ^ S (6 + S) ^ 5 (S - if 
Tl 1 + 6 S 7 1 + 6 S 


= Ms - 1) 

Tl 7 


(^) 


(5) 


where S > 1. Relation ( 5 ) gives the theoretical rise in stagnation 
temx>erature behind the primary shock wave traveling in still air. The 
plot of Tq 2 versus S (also for other temperatures of possible interest) 

is given in figure 5 * 
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Contact -Surface Region (state c) 

State c Is the region that contains the gas ^ich^ to start with, 
was in the compression chamber. As figure 2 shows, it ends at the foot 
of the centered rarefaction wave. The contact-surface flow is of con- 
siderable interest in shock-tube work. The sudden expansion of the gas 
in the compression chamber cools the gas considerably, thus reducing the 
local speed of sound, which in turn yields higher Mach number flows. 
Assuming that there is no change in entropy through the rarefaction wave 
generated because of this sudden expansion, one can solve the remaining 
continuity and momentum equations by the method of characteristics 
(ref. 11) such that the unknown physical quantities through the wave can 
be related to the known quantities ahead of the wave. If the centered 
rarefaction wave moving to the left in a gas at rest (shock-tube case) 
is considered under the assumptions that the expansion of the gas in the 
ccxnpression chamber is adiabatic and 


Pc = P2 

> 

Uc = u^^ 

then for Air /Air operation. 


P = ^ = s 
Pi 


1 


S - 1 

\Zt(i + 6s )_ 


-7 


( 6 ) 


(7) 


The derivation of this basic equation is given in almost all of the shock- 
tube references listed. The physical quantities such as temperature, 
press\ire, density, speed of sound, wave speeds, and particle velocities 
throughout regions r, c, and 2 are related to known quantities in 
regions R and 1 through relation (j) which connects the pressxire ratio 
across the diaphragm to the shock strength. From the consideration of 

Mach nvnnber = — , and making use of the boiindary conditions (6) and 

&c 

the fact the expansion is assumed adiabatic, one can get 


Tr 


‘■R 


-.2 


S - 1 


^(a + 1)(1 + aS) 


( 8 ) 


It has been experimentally observed that the hot-wire responds to the 
stagnation temperature of the flow. With that end in view, an expression 
for the stagnation temperature in the contact -surface flow can be obtained 


12 


NAGA TN 5163 


from the consideration of the energy equation of the contact -surface 
flow. One gets 


*^00 _ *^00 _ 
Ti 'Tr 


1 + 




- 1 


KR 


- 1 


1 - 


S - 1 


- .ih r. - 


\J(c + l)(l + crS) 


S - L 


^(a + l)(l + aS) 


(9) 


Plots Of Tf. and Tq,, versus S are given in figirre 5. 


In calculating hot-wire voltage signals for piore temperature response, 
one has to consider the fact that the hot-wire reaches equilibrium at 
approximately the stagnation temperat\ire of the flow behind the primary 
shock wave before the contact surface washes over it. So the decrease in 
resistance due to the cooling of the wire will start from the elevated 
temperature. Therefore Tq 2 - Tq^, is more useful. From relations (5) 

and ( 9 ) one can get 



S - I y 

5 _ - 

i (s - 1) + 1 

^/ 7(1 + 6S) 

7 


1 


(10) 


One also can obtain the mass flow in the contact surface given by 


Pc^c _ p 

1 - s - 1 , . 1 

(a - 1)(S - 1) 

Pl^l 

]j{a + 1)(1 + aS) 

^(a + l)(l + oS) 


( 11 ) 


A plot of this is given in figure 4. For an Air/Air operation the ratio 
of relation (ll) to (2) is pc/p2- consideration of the ordinary 

gas equation 


^ = 


P2 



Since T 2 > Pc ^ ^2" there is no mass flow across the 

contact -surface discontinuity, it involves a sudden temperature and den- 
sity discontinuity between the mass of the gas which originally belonged 
on the two sides of the diaphragm. This fact has an Interesting bearing 
on the Reynolds numbers in adjacent regions 2 and c of the flow (figs. 6 
and 7 ) • 
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Reflected Shock and Region Behind It (state 3) 

In the conventional arrangement in shock-tuhe vork^ the primaiy shock 
wave travels through the medium which is at rest in front of it. If the 
far end of the shock tube is closed by a completely solid and rigid wall 
parallel to the shock front, the shock will be totally reflected and will 
travel upstream against the drift velocity of the primary (incident) shock 
wave. The gas in contact with the rigid wall must be at rest. The 
reflected shock thus has a strength such that the drift velocity due to 
the primary shock heading toward the closed end is canceled by the reflec- 
ted shock. While the primary shock accelerates the medium from the speed 
zero to U 2 ^ the reflected shock decelerates it back to zero. Thus the 
reflected shock differs from the primary shock in that the state in front 
of the reflected shock wave is identical with the state behind the pri- 
mary shock wave. Let the reflected shock strength S' = PWP 2 - From the 
Rankine-Hugoniot relations for the incident and reflected snock and the 
condition that 


U2 = U2* 


(where primed parameters are associated with the reflected shock) one can 
obtain the speed of the reflected shock with respect to the shock tube as 


U 


s * 


ai 


2S + (g - 1) 
\/(g + l)(gS + 1) 


The temperature behind the reflected shock is 

^ = S(o + S) (g + 2)S - 1 / 2 ^ g - 1 l\ 

1 + gS g + S \g + 1 g + 1 S/ 

A plot of this is given in figure 5* 


( 12 ) 


Rarefaction Waves in Shock Tube 

Centered rarefaction wave propagating in compression chamber 
(region r) .- From the viewpoint of the hot-wire anemometer as a timing 
device, the centered rarefaction wave (region r) is also of some interest. 
It originates because of the sudden adiabatic expansion at the abrupt 
breaking of the diaphragm between the expansion and compression chamber. 
Such a wave represents a continuous disturbance for t > 0 where t = 0 
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is the time of rupture. If the problem is analyzed by the method of 
characteristics^ one can arrive at the following conclusions: 

( 1 ) The head of the rarefaction wave travels with the speed of sound 
of the medium at rest in the compression chamber. 

( 2 ) The temperature change across the rarefaction front and the mass 
flow per unit area initiated by a centered rarefaction are given by 

T = ^ 

and 

Pr^ + 1 V'T /\7r + 1 °R 

where 


(15) 


(i4) 



T 


fR^ 

L 


and L is the length of the compression chamber. For the desired hot- 
wire use as a timing device, *^r/^ Pr^r/PR^R main interest. 

These are plotted in figures 8 and 9 for various times after the diaphragm 
is ruptured, namely, for t = 0.20, 0.25, O. 5 , O. 75 , and 1.0, against 
various locations in the shock tube. Interaction of these plots with the 
corresponding temperature and mass flows in the contact surface (rela- 
tions ( 8 ) and ( 11 ), respectively) sets the limit to 
application of relations ( 13 ) and (l4). 

For the head of the rarefaction wave, at = 


the region of 


X 


a^t 


= -1, one gets 




Similarly, 


(15) 


^ 2 

+ 1 


(16) 
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Thp head of the rarefaction thus represents a discontinuity in tempera- 
ture and mass flow. The hot-wire responds to both temperature and mass- 
flow change. The rsLref action initiates cooler flow. If the hot-wire 
is used at an elevated temperatiire, both the mass-flow and temperature 
changes involved give signals of the same sign. Thus one can use the 
hot-wire technique to clock the speed of the rarefaction wave. Since 
the head of the rarefaction travels at the speed of sound in the still 
gas in the compression chamber, one can measure the sound speeds in 
various gases rather easily with the hot-wire technique provided the 
hot-wire is used at an elevated temperature. 

Rarefaction originating at open end of shock tube .- If the shock 
tube is operated with the extreme end of the expansion chamber open 
(fig. 2), a rarefaction wave originates when the primary shock wave 
leaves the open end. This is due to the pressure behind the traveling 
primary shock wave being higher than the atmospheric pressure surrounding 
the shock tube. If the flow behind the primary shock wave is supersonic, 
the open end does not influence the quasi-steady flow inside the shock 
tube. But for the subsonic range (l.e., M2<l)a rarefaction wave is 
propagated toward the inside of the shock tube starting at the instant 
the primary shock wave leaves the open end (fig. 2 ). This wave is, in 
general, transmitted as a rarefaction wave, at least until it encoxmters 
the contact siarface inside the shock tube. The propagation of this rare- 
faction wave has been studied experimentally (procedure and results to 
be discussed later), especially from the point of view of its speed, that 
is, the time of its appearance at fixed stations along the axis of the 
shock tube under different working pressure ratios, and so forth. For 
the detection of the time of arrival, it is sufficient to consider only 
the first characteristic of the wave originating at the open end. It 
will fiorther be assumed that the flow behind the primary shock is sub- 
sonic that the position of the detecting station is at such a distance 
from the open end that the contact-surface flow or the principal reflected 
rarefaction wave does not reach that station prior to this rarefaction 
front under investigation. The first characteristic of this rarefaction 
wave travels at the speed of sound in the region of transmission. Based 
on these assumptions the time of arrival of the first characteristic at 
any station can be predicted. 

Suppose I and t] are the coordinates of the detecting station 
(which, in this case, will be a hot-wire moimted along the axis of the 
shock tube). Then, as shown in figure 10 , 



(17) 
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y = n 


(18) 


X = I + u^t 


where U2 is the particle speed and a2 is the speed of soiond in 
region 2 . 

Therefore from relations (17) and (18), the velocity of propagation 
of the rarefaction front moving towaxd the inside from the open end of 
the shock tuhe is 


Retaining the sign which gives positive velocity of propagation, for 

0 = 0 , 


For the outside of the shock tube, a similar procediire gives the velocity 
of propagation as 


Retaining the sign which gives the positive velocity, for 0 = 0 , 




( 19 ) 



^ri - - U2 


(20) 



( 21 ) 



( 22 ) 


For 0 = 90°^ that is, at the open end, frcm relations (l 9 ) and ( 21 ) 


Vj. = ^ sl 2 ^ - ^ 2 ^ 


(23) 
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If a hot-wire m.o\inted in the middle of the shock tube as a detecting 
device is used, then, when the primary shock wave passes over the hot- 
wire, it will register the jump and the response will stay steady until 
the rarefaction wave passes over it (fig. l). The total time interval 
between the instant of jump due to the shock wave and the time the 
cooling due to rarefaction is registered will be 




u. 


Vri 


(24) 


If the hot-wire was mounted in the wall of the shock tube, then 


= -L + _i_ 
Us V^i 


(25) 


For the hot-wire mounted outside the shock tube along the axis 


t = 


J_ 

Vro ‘ Us 


(26) 


where | is the horizontal distance from the open end. (Experimental 
verification of these time intervals will be discussed later.) 

The rarefaction wave completely overtakes the whole shock front at 
about an equivalent length of a diameter of the shock tube downstream 
from the open end. The decay of the shock front becomes very evident 
after that. 


SCME RELEVANT FEATURES OF HOT-WIRE RESPONSE 
Steady-State Response of a Hot-Wire 


If a wire, through which some heating current is passing, is 
immersed in a flowing fluid, the wire responds to the changes in the 
heating current or in the cooling (or heating) effects of the flow 
according to certain response equations. 

In order to derive these equations the conservation of energy is 
assumed, that is, the difference between electric heat input and heat 
loss due to the forced convection (including radiation and free convec- 
tion, etc.) eqtials the increased thermal energy in the wire: 
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W 


H = C 


dT, 

dt 


V 


where 

H heat loss from wire per unit time 

W energy input in the wire per unit time 

For the steady case 


(27) 


W - H = 0 


(28) 


The well-known investigation of heat loss from the cylindrical wire placed 
in an incompressible (M — > 0) flow conducted by King (ref. 12) gives for 
heat loss 

H = H'l’w - Te) (ke + ^2«kgCppud) 


or 


H = (T^ - Tg)(A + B^) 


( 29 ) 


where 


A = kei 
B — i ^ 2 itk^Cppd 


Using relation ( 28 ) 

I^Rw = (Tw - Te)(A + B(H) 


This can he further transformed into (refs. 4, 5^ sind 10) 

lo^ ^ 


(50) 


( 51 ) 
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This is a steady-state equation which governs the hot-wire heat-loss 
behavior in a steady stream of flow. The characteristic consteints Iq 

and Uq are functions of temperature and pressure or density of the 
flow and are experimentally determined by the calibration of the wire 
(see section "Calibrating Facility"). If the flow conditions of the 
calibrating jet are different from those of the flow under investigation, 
then Iq and Uq either will need modification or, in more extreme 
cases, may exhaust their usefulness and thus may have to be discarded 
for an alternative approach. Qlie variation of these constants is usually 
smal 1 if the operating conditions are not radically different from the 
calibrating conditions. The variation of Uq due to the actual flow 
being of a substantially different density cam be taken into account 
rather simply by expressing relation ( 51 ) in the form 



(52) 


where is some reference ^ density and p is the density of the 

actual flow. Supposing Uq = 10 feet i>er second at normal atmospheric 

density, then, if the density of the flow was, say, one-half the normal 
atmospheric density, Uq will be 20 feet per second. The function Uq 
does not depend upon the material of the wire and for constant tempera- 
ture Uod = Constant. The characteristic constant Iq is different for 
different types of wires (fig. 11). Since in shock-tube flows the temper- 
atures vary considerably from the normal room temperature, the dependence 
of k^ on temperature has to be taken into account. (Variations of Cp 
with moderately high temperatures may be neglected.) It is sometimes 
advantageous to put the heat-loss equation ( 51 ) into an alternative form 
valid under widely varying conditions: 


Nu = a + b^|p^ 


(55) 


where a and b are constants and the Peclet number Pe = PrRe. When 
centimeter-gram-second units are used then Nusselt's number 

^ 0.2389 I^Ry, 

rtZkg^Tw - Te^ 
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Therefore, 


0.2389 I2Rv 


= 0.318 + 0.798/^ 


jt Zkg 


If Pr = 0.72, which for diatomic gases is a close value over a wide 
range of temperatures (ref. 13 )> then 


Besides the forced convection, the heat loss is influenced quite appre- 
ciably by free convection, radiation, and conduction thro\igh ends. 

For hot-wires used in low-density wind tunnels or low-density opera- 
tion of shock tubes, the heat loss due to free convection is substantial. 
If the hot-wire is heated to a certain temperature under normal atmos- 
pheric conditions, then much less current is needed to heat the same wire 
to the same temperature \mder partial vacu\jm conditions. To calibrate 
this effect a platinum wire (d = 0.0001 inch) and a tungsten wire 
(d = 0.00015 inch) were mounted in a shock tube. They were heated to 
a^ = 0.5 and a^ = 0.1 at various pressures within the range 760 
to 0.15 millimeter of mercury. The corresponding cold resistance was 
measured at each pressure and the operating resistance calculated for 
the different overheating ratios. The ciorrents needed to heat the .wires 
to these overheating ratios were measvired. These wires were mounted both 
vertical and horizontal to the base of the shock tube . Since no appre- 
ciable difference was noticed between the performance of these two orien- 
tations of the wires, only the data for vertical wires are reproduced 
here . 

Even in still air, there is always heat loss. This may be due to 
free convection or molecular conduction. 


The measure of free convection is the so-called Grashof number 
( ref . l 4 ) 


where g is the acceleration due to gravity and Tg is the gas tempera- 
ture surrounding the wire. 



( 54 ) 


Gr = 
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Frcm relation (55) (for 
the gas 

- T, 


This could he written as 


%°-l j2 
1 


If the temperature of the gas 

j2 1+ a-M ^ a(awP^) (57) 

The examination of the data from the point of view that free convection 
is entirely influenced by Grashof number failed to correlate the data 
for a^ = 0.5 and = 0.1 (see ref. 10 ). Because of the smallness 
of the diameter of the wire, the Grashof number even at atmospheric pres 
svire is very small in magnitude (approximately 10“°) and the normal free 
convection effect (with usual free -convection velocity, and the associ- 
ated buoyancy effects) is considerably reduced (ref. l4). 

Free convection is also affected by the moleciilar conductivity 
through the dimensionless quantity known as the Knudsen mmiber 


= 0 ) and for constant temperature of 


oc a 


gp^d.^(Tw “ Te) 




(55) 


1 + ^w 
aw 


<x a 


'gp^d^a^ 
M- o^e^e 


( 56 ) 


is constant, then for the same wire 



where 

A molecular mean free path 

d characteristic length (diameter of wire) 

Since A is inversely proportional to the density of the gas, for the 
same wire, K is greater at lower densities. This shows no influence 
of overheating ratio hut does involve density. 
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Relation (35) can "be written as 

^ fxjLl ccaW (58) 

a-y \d/ 

where a is the specific resistivity of the material of the wire. The 
data are plotted as shown in figure 12 and the values of d, oe, o, 
and A have been assumed as Indicated. 

The decrease in the amount of current needed to heat the wire to the 
same temperature under partial vacuum is striking. 

The heat-loss behavior of both tvingsten and platinum wires agrees 
for large values of K (low densities), that is, where the mean free 
path is large compared with the characteristic length. 

For comparatively higher densities, however, the agreement for 
platinum wire for the two overheating ratios is poor. This disagreement 
is not fully \mderstood. 


Linearized Unsteady Response to Fluctuations 


The approach used is that of reference 10. From the consideration 
of energy conservation, relation (27) was written. If one assumes that 
the flow cam be split into mean motion and the fluctuations (A quanti- 
ties), and since mean flow satisfies the steady-state relation (28), 
one gets 


AW - AH = 


p <i(^w) 
dt 


( 59 ) 


where 


ah = AT. 


8t 


W 


' * ' a(pu) 


A(pu) 
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Collecting and dividing by the coefficients of and from the knowl- 

edge of the explicit expressions for W and H, one can obtain the 
expression 


M* = 

^w ' ^w 



i^O) 


where 


n 


C 

oti^ 


aw = 


Rw - 
Re 


The thermal-lag constant n is a universal constant of the wire as it 
is directly related to its physical properties. For a specific material 
it varies with the f ovirth power of the diameter of the wire . 

The voltage fluct\aations across the wire terminals can be due to 
resistance and/or current fluctiiations 


Ae = I AR + R^ AI 


( 1 + 1 ) 


Using relation (1+1) along with the assumed heat-loss law (King's incom- 
pressible case)^ relation (39) can be written as 


-F A(pu) + G ATe 


^ + M** 

dt 


(1+2) 


where F and G, the sensitivity coefficients, and M** are given by 
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TO _ eaiEZ 

G = oEjIE ^ > 


(43) 


M** 


M* 

1 + 2a^€ 




and 



Here 


(pu)o is the calibration constant and 


= 1 - c 

1 + 2a^e 


where 


e = 




+ Rw 


"open 


Zg impedance of current source 

e mean voltage drop across the wire, 1^ 

e open-circuit voltage of hot-wire supply circuit if hot-wire 

open ^ • 

is removed 

If the alternating-current source impedance of the hot-wire terminals 
is infinite, then e — ^ 0 and E = 1 and therefore 
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F = ea^z/k!pu 


M* = na^/l^ 


m 


G — ZiR^HxR-^ 


J 


If IV 


Rg , then 


G = loRf 


(^5) 


which is the case when the current passed throiigh the wire placed in the 
flow does not raise the temperat\ire or resistance of the wire appreciably 
higher than the equilibrium temperature which the wire assumes without 
any current through it. 

If the conditions vary extremely slowly or the wire has no themal 
heat capacity, then the left-hand side of equation (42) is the static 
response of the wire and one csua write 


where is the virtual voltage fluctuation which the wire, without 

any thermal lag, would have assumed and M* is the usual time constant 
of the wire at particular operating conditions. Equation (46) is equiv- 
alent to the differential equation obeyed by a circuit containing an 
inductance and resistance or resistance and capacitance arranged as an 
integrating device. Such equivalent circuits can be used as "dummy hot- 
wires" (ref. 4). 


The principal application of hot-wires is to turbulence measurements 
where the fluctuations usually are small compared with the mean motion 
and the unsteady hot-wire response equation is linearized resulting in 
a time constant M* which is a function of the average operating con- 
ditions (relation (44)). If a traveling shock or any other transient 
wave passes over a hot-wire placed in the shock tube, both the mass flow 
and the ambient temperature change as a step function in time. This fact 


= ^ + M* 


d(/!^) 

dt 


(h6) 


Response to a Finite Step Function 
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enables one to solve the exact equation without resorting to the pertur- 
bation method and results in a simple rule that the time constant may be 
determined by the operating conditions that would be assumed in a steady 
state a^ter the jvunp. 

For this type of derivation the following assumptions are made: 

(1) The current I remains constant (reasonably well justified by 
the hot-wire equipment used). 

(2) Thermal condition of the wire can be characterized by a single 
temperature T^ (thermal wave propagation along the wire is negligible). 
This means that the distribution shape of the temperature along the wire 
does not change (except for a constant) ixnder different operating 
conditions . 

(3) The shock fronts, contact surface, and reflected shock front are 
all plane. 

( 4 ) King's heat-loss law applies for the incompressible range (low 
Mach numbers) and Kovasznay's empirical heat-loss relation applies for 
the supersonic range considered. (The limitations of this assumption 
will be dealt with later.) 

(5) The temperature which the wire will assume if it was not heated 
at all will be called the equilibrium temperature and it is approximately 
equal to stagnation temperature. It is assumed that the hot-wire responds 
to stagnation temperature. 

Primary shock wave.- The heat generated at any time is 



(^ 7 ) 


The heat loss after the primary shock has passed is 



m 


The hot-wire transient response equation (27) becomes 



(Tw - To2)(a + = C ^ 


()+9) 


NACA TN 5163 


27 


After equilibrium between wire and flow has been established. 



and Ty— 


0 




1 + 


<‘^2 - Tf) - (Tv2 - = ° (50) 


Subtracting equation (50) from {kS) and eliminating A + one 

gets 


T^2 - T. + m/ ^ (51) 

where 

T _ T 

M * = - 9 — 

^ I%f 1 + a(To2 - Tf) 

= ^ a ^2 (52) 


with 


%2 = 


^2 ~ ^o2 


R 


o2 


In this special case when the mass flow P2U2 and temperature Tq 2 
constant in time, except for discontinuous jumps from the ambient con- 
dition, M2* still remains constant except that its value must be com- 
puted with the steady-state operating conditions which the wire assumes 
"long" after the transient occvirs. This means that wire response will 
still be an exponential function of time and that ordinary electronic 
thermal-lag compensation sho\ild give proper compensation if the time 


constant is set according to 



This has been checked to be 


so as described in reference 8 and rechecked again (see "Analysis of 
Experimental Data"). A typical compensated and uncompensated oscillo- 
gram is reproduced in figure 1 . If the fluctuations in the flow behind 
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this main jiomp are being studied, this relation for compensations is 
applicable only if the fluctuations are small as compared with the mean 
flow. 


Contact surface.- After the shock has passed, the wire adjusts to 
the conditions behind the shock and when the contact discontinuity reaches 
the hot-wire it experiences another step-function change in flow condi- 
tions. The temperature of the wire before the contact surface reaches it 
is Ty2> and after equilibrium with the contact surface is established 
the wire temperature becomes Following a procedure similar to the 

one used for the derivation of equation (^ 1 ), one gets 


VC 


= T o + M * 
•^v2 c 


dT. 


v2 


dt 


(55) 


where 


Me* = 



( 5 ^) 


with 




R. 


VC 


- R, 


oc 


R, 


oc 


Comparing relations (52) and (5^)^ 

(55) 

^c* ^•vc 

where ay 2 and aye are the new equilibriiam overheating ratios which 
the wire acquires in the flow behind shock and contact surface, respec- 
tively. From the calculations (details to follow in the next section), 
8 y 2 > aye by approximately 20 to 30 percent for the same shock 

strength S. This means that lesser compensations will be needed for 
the contact -surface region. 
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So far the discussion has centered around the Jumps due to the main 
wave fronts (shock waves, contact surface, and rarefaction wave). If 
there are fluctuations present (natural or artificially introduced) in 
the flows behind the main fronts, the linearized theory will hold again, 
with the difference that the new equilibrium operating conditions a^2 

and a^^ are different from the original a^. This affects the sensi- 
tivity coefficients of the hot-wire response (relations (kk) and ( 45 )). 

COMPUTATION OF EXPECTED VOLTAGE SIGNALS FROM HOT-WIRE 
OPERATED IN SHOCK TUBES 


The contents of this section are in general helpful for the opera- 
tion of the hot-wire anemometer and the design of the associated equip- 
ment for use in the shock- tube work. It gives the order of magnitude of 
voltage signals expected from the various regions of the flow under some 
standard operating conditions. Until otherwise stated, the following 
assumptions have been made: 

(1) For the subsonic range of flow, King*s heat-loss relation is 
considered valid. This means that Mach number effects on the heat-loss 
behavior of the hot-wire used in subsonic flows are overlooked. 

(2) The wire responds to stagnation temperatures. 

(5) The hot-wire is heated from a current source of infinite 
impedance . 


Region Behind Primary Shock 

The flow behind a shock wave has temperature and density higher than 
that of the region ahead of it. Since the calibration constants of the 
hot-wire Uq or (pu)o and 1 ^ are functions of temperature, for high- 
strength shocks they have to be modified. The alternative, of course, is 
to resort to the actual physical dimensions and so forth of the wire; 
while for this case one does not need to calibrate the wire, one has to 
rely on the measured length and diameter of the wire which are rather 
poorly defined and vary from sample to sample. So to calculate the jump 
signals due to response of the hot-wire to a traveling shock wave, one 
is offered a relative choice at best. Calculations carried out by both 
methods did not show any pronounced difference. Since the aim of these 
calculations is only the estimate of the order of magnitude of the sig- 
nals involved, the choice is not critical. The approach elaborated here 
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is the one using the physical dimensions of the wire. In conventional 
shock-tube work there is no flow to start with, that is, \f^ = 0. 
Therefore, from relation (54), the heat loss for no flow is given in 
centimeter-gram-second units by 


Nu^ 


0.2389 

nlke^Ty - Tgj 


0.318 


(56) 


Since because of the passage of the shock wave \f^ jumps to yR^, 
one gets 


2 

^ — 0-2589 1 .... % — ^ ^ 0.318 + 0.6771/Re^ (57) 

itZko2(Tv2 - T 02 ) 


where 


Re2 = P2^2^/^^o2 


For constant -current operation 


and 


Nu2 ^v 2 *^w " 


^ ^o2 *^v2 *^o2 


where 


_ ^v2 ^w ~ 

^2 p ^ 


^w ^v2 “ 


o2 


Po - 


Nu2 ^q2 
Nu 3_ 


therefore 


Rv 2 = 


1 


1^02 

1 ^ 

P2 &w + 1 


(58) 


(59) 
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From relation (56), 


^2 0.5l8jfk.]^^ay 

“ 0.2389RwO'1 


(60) 


If the respective values of Z and are selected, then for various 

operating overheating ratios, the corresponding currents can be calcu- 
lated. The values assumed for the tiuigsten wire of d = 0.00015 inch 
are = 10 ohms, I = 0.2 centimeter, and 

= 5.71 X 10“^ cal cm"^ sec"^ °C"^. If the shock strength S is knoi 
then ^Re2 and Nu2 can be calciolated and, since Tq 2 is known from 
relation (5), corresponding air conductivities kQ2 are obtained from 
reference 15 and P2 calculated from relation (58). Since Tq2 

is known, Rq 2 can be found. As a^ is known at the start, Ry2 can 
be calciilated from relation (59)- Also ay2 can be calculated. The 
voltage jump for constant-ciirrent operation is 


where tuB is positive if R^ < Ry2 ®cid negative if R^ > Rv2* 

The shock tube can be operated in many different ways. For instance, 
one can keep the expansion chamber at normal atmospheric pressure and 
pressurize the compression chamber or one may partially evacuate the 
expansion chamber and keep the compression chamber at normal atmospheric 
conditions (or higher). The method of operation has a striking effect on 
the \/Re of the flows (figs. 6 and 7). 

The calculations of Ae's as a function of shock strengths with a^ 
as a parameter presented here assume that the shock tube is operated with 
the expansion chamber at normal atmospheric conditions. Notice that there 
are combinations of a^ and S where no signal is available (fig. 13)- 
The shock wave represents a mass-flow as well as a temperature jump and 
the hot-wire is sensitive to both these changes. While the temperature 
j ump tends to increase the resistance, the mass flow reduces it. These 
signals are of opposite polarity. Thus in between there must be a 
"neutral" point where these two signals cancel out, giving no response 
at all. The point is well borne out by the oscillogram produced in 
figure Ik. Corresponding to the Intermediate cvirrent there is hardly any 
jump, while the largest and the smallest currents have jumps of the oppo- 
site polarity. Needless to say, while the hot-wire is used for timing 



(61) 
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p\irposes, one shovild avoid working in the neighborhood of this region. 

It is advantageous, as a matter of fact, to stay as far away from it as 
possible. However, when the wire is used for actual flow studies this 
presents very little difficulty. On the contrary, at times it may be 
advantageous to operate closer to this point in order to avoid large 
voltage signals. Notice the extremely large signals obtainable by oper- 
ating the hot-wire with a considerable current throxigh it at extremely 
weak shock waves. In the author's opinion this is the region where the 
hot-wire is uniquely we3_L sxxlted for timing shock-wave speeds, and so 
forth. 

During the course of the work described in this report by chance a 
tungsten wire was found with a resistance of exactly 10 ohms. The wire 
was calibrated and actioally used for measuring jump voltages for low- 
strength shocks at overheating ratios of 0.5 and 0.3- The expansion 
chamber was at atmospheric pressure. The amoimt of signals obtained 
(fig. 13 (a)) verifies the order-of -magnitude ceilculations for the jump 
signals . 

Tungsten wires oxidize at much lower temperatures than platinum 
wires. Based on personal experience, the author would advise against 
the use of t\ingsten wires beyond an overheating ratio = 0.7- For 
very weak shocks platinum wires are more desirable since one can beat 
them to higher overheating ratios without their oxidizing. 


Contact-Surface Flow 

The contact-surfact flow follows the flow behind the primary shock 
front, thus introducing another sudden change in temperature, mass flow, 
and Mach number- The hot-wire comes into equilibrium with this flow. 

But in this case the ini tial conditions are the equilibrixmi conditions 
due to the flow behind the primary shoc k wave . The Reynolds numbers of 
the hot-wire jump frcan )/Re 2 to ifRec- As the Mach nimiber effects 
(Me < 1) wi3J- again be neglected, the corresponding Nu^ is calculated 
from King's heat-loss eqioation. By a similar procedure to that followed 
for the flow behind the shock wave. 



Qy2 

Pc 1 + ®v2 


(62) 




where 
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Knowing and Tqc can find Rqc sind thus Rye can be deter- 

mined from relation ( 62 ), and the voltage Jimp is given by 


= I^Rv2 " ^c) 


(63) 


Since Ry2 > Rye these signals are negative and increase with the 
increasing overheating ratio initially used (fig. 13(t>)). The values 

of used here are the ones calculated with the expansion chamber 

at normal atmospheric conditions (ref. 9). 

If one wants to use the contact surface as the principal front to 
start the triggering circuit, one should operate it such that the shock- 
wave signal is positive in polarity and the associated electronic equip- 
ment should be so arranged that it overlooks the signal due to the shock 
wave and does not trip the thyratron tube until the negative signal from 
the contact svirface activates the system (ref. 9)* 


Pure Temperature Response 

From the consideration of the sensitivity of the hot-wire (rela- 
tions {hh) and (4-5)) one observes that, if the operating resistance of 
the hot-wire is very close to the equilibrium resistance of the wire 
without any current through it, the temperat\ire sensitivity dominates 
and the hot-wire acts like a resistance thermometer. If the operating 
temi)erattire of the wire is not much different from the gas surroxmdlng 
it (i.e., "small" current) and if the current is assumed constsmt, then 
the voltage change is given by 


Ae = ±I AR 

= tlRjLa]^ AT 


(64) 


where AT is the change in temperature of the gas passing over the hot- 
wire. If the resistance of the wire is, say, 10 ohms, the current through 
it is 4 milllamperes, and a = 0.0035 then the voltage Jump per degree 
Jump in temperature is Ae = 0.l4 mv/°C. The temperature Jumps across 
the shock wave, contact surface, and reflected shock are of such magni- 
tudes (fig. 5 ) that one can obtain fairly large voltage signals from the 
timing hot-wire even for moderately strong shock waves. 

However, the temperature Jumps across the strong waves are so large 
that, even if one operates the wire at relatively high temperatures, the 
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temperature response dominates and the signal is positive. The oscillo- 
graphic record of the hot-wire response reproduced in figure 15 demon- 
strates this clearly. The initial jump (A) is due to the shock wave and 
after a millisecond or so the contact surface arrives at the hot-wire 
and cools it. The tiirbulent nature of the contact-surface flow (b) is 
evident . 

The break in the first Jump is due to the firing of the spark (D), 
the output of which was pxirposely applied to the oscilloscope. 


HOT-WIRE RESPONSE TO TRANSIENT SUPERSONIC FLOW 
REGIONS IN SHOCK TUBE 


In shock-tube investigations, the transient supersonic flows are 
very easily produced. There are two regions of such flows which are of 
main interest: 

(1) Flow behind the primary shock wave (region 2 ). 

(2) Flow in the contact surface (region c). 

While the pressure and particle velocity are the same in both the regions, 
the temperature and density are quite different. Therefore, for the same 
shock strength the Reynolds numbers suad Mach numbers in these regions 
differ quite substantially. 

Hot-wire techniques have been recently extended to supersonic flows. 
Kovasznay’s empirical law of heat loss obeyed by the hot-wire placed in 
a supersonic stream is given by 



where A, B, and C are nondimens ional constants. For low- temperature 
loading 
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The heat loss is still dependent on the Reynolds mamber. The negative 
intercept in the plot of Nu versus (ref. I6) is different from 

the low-speed heat-loss la.w (King's case). Experimental stiodies of heat 
transfer from bodies in a high-speed rarefied gas stream (i.e., Reynolds 
number lower than the rarige at which Kovasznay worked) show that this 
negative intercept in reality does not occiir. Starting near the origin 

Nu versus ^e is parabolic in this lower range of Re and joins with 
Kovasznay 's empirical curve (l.e., if fig. 7 iii ref. I7 is replotted as 
Nu versus ) . 

For the study of heat transfer from bodies in high-speed and low- 
density flows, in addition to Reynolds number and Mach number, the 
IQiudsen number K = A/d becomes Importemt, where A is the molecular 
mean free path and d is the characteristic length (diameter of the 
wire in the present case). 

The parameter K can be related to M and Re as (refs. I7 and 

18) 


K 



One usually uses the Knudsen number as a criterion for deciding in which 
type of flow one is operating (l.e., continuum, slip, or free moleciolar). 
When the Knudsen number is small (K < 0.01) the effect of the molecular 
motion on the flow is negligible and the gas may be treated as a continuum. 
For large values of Knudsen mmiber (K ^ 10) the moleciilar motion is all 
Important and the phenomena can be completely described in terms of free 
molecules. The realm of flow with Knudsen mmiber 0.01 < K < 10 is 
usually classified as the slip-flow region. In such a flow the gas no 
longer sticks to the surface but slips over the surface with definite 
velocity. 

The shock tube can be operated in a number of ways. The method of 
operation has rather a striking effect on the Reynolds number obtained. 

This in turn has an effect on the Knudsen number and the heat loss from 
the wire. Calculations of Reynolds numbers and Knudsen numbers were 
carried through for the flow behind the shock wave and contact surface 
for the following cases: 

(1) The expansion chamber at atmospheric conditions and the compres- 
sion chamber at higher pressures (figs. 6 and 16). 

(2) The expansion chamber partially evacuated and the compression 
chamber at atmospheric pressure (figs. 7 and I7). 
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A few salient features are stated here . 




He c • 

(a) The Reynolds ntunber in the contact siirface Re^, = is 

M-oc 

higher than the corresponding one in the region behind the shock 

^ This is due to the higher density and lower temperature 


^^02 


Re2 = 

(i.e., lower viscosity coefficient ref. 19) fJi the contact-surface 

flow. 


(b) When the shock tube is operated with the expansion chamber open 
to the atmosphere, the Reynolds numbers are higher than the corresponding 
ones for the expansion chamber operated at low pressures. 

(c) Based on the criterion discussed above, the numerical values of 
the Knudsen number show that it is only when the shock tube is operated 
with the expansion chamber at normal atmospheric pressure that the flow 
over the hot-wire can be treated as a continuum. When the expansion 
chamber is operated at low densities, K is small and one operates the 
jjot— wire in or very close to the slip— flow conditions for a large range 
of shock strength. 

Because of the small characteristic length (diameter of the wire 
d = 0.00015 inch), even at normal temperature and pressure, K = 0.02 

-6 

with A = 8.5 X 10" centimeter. 

Most of the turbulence investigations conducted with the hot-wire 
technique in the main stream of a supersonic wind tunnel have 
K « 0.02 — ^ 0 . 0 k. This is also the range in which Kovasznay worked to 
establish the empirical heat-loss law (relation (65)) which breaks down 
for K 0 . 1 . According to the cirteria agreed upon, this lies in the 
slip-flow region. 

Thus while operating the hot-wire in shock tubes one should decide 
the heat-loss relation to be used (and therefore the estimated signals, 
etc.) according to the operating range of Reynolds number, Mach number, 
and Knudsen number. 


EXPERIMENTAL EQUIPMENT AND TECHNIQUE 


The experimental arrangement consists of the following f\mctional 
units: 

(1) Shock tube 

(2) Hot-wire equipment 
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(3) Triggering and timing 

(U) Mounting of hot-wires 

( 5 ) Recording 

(6) Calibrating facility 

Relevant featvires of each unit are summarized in the subsections which 
foiJLow. 

Shock Tube 

The shock tube consists of a steel box of \mlform cross section 
(approximately 6 by 4 inches inside) with the 6 -inch side laid horizontal 
(fig. 18 ) to allow the largest path of the beam of light throxigh the flow 
\mder Investigation. The total length of the shock tube is 29 feet. The 
ccanpression chamber is 6 feet long and the expansion chamber consists of 
seven sections (fig* 19)* 

The construction of the shock tube is weld free and, as this was 
intended to be a functional rather than a precision shock tube, the 
minimnm amount of machining was employed. Most of the shock tube is 
constructed of ordinary cold-rolled steel plates. The plates are bolted 
together with l/l 6 -inch rubber gasket as sealing material. Dowel pins 
are used to keep the box square and alined (fig. I 8 ). 

The construction of the glass section is somewhat different. This 
is the only completely machined part in the whole shock tube. The glass 
window extends from the top to the bottom plate and one can have any 
section of the glass free from bolts with a maximum clear view of 4 by 
11 inches. 

The ccanpression chamber is capable of withstanding considerable 
internal pressures. It has been used at 125 pounds per square inch. 

The high pressure is measured with a 12-inch-diameter pressure gage 
gradviated from 0 to 25 O pounds per square inch. 

When the expansion chamber is used \mder vacuum the pressure is mea- 
sured by a merciiry manometer, absolute pressure gages (O to 20 millimeters 
0 to 100 millimeters, Wallace and Tieman), and a McLeod gage (O.Ol to 
50 millimeters). Among them the full range of low-pressure measurement, 
which one may feasibly require in shock-tube work, is covered. The over- 
lapping of their ranges is used to check one against the other. 

The shock tube is provided with a spring-loaded, solenoid-activated 
mechanism with a thin pointed needle to punctiare the diaphragm. 
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General performance of this simply constructed shock tube is very 
satisfactory whether it is used at high pressures or with expansion under 
partial vacuum. The expansion chamber has been used at pressures lower 
than 1/10 millimeter of mercury. (For details see ref. 10.) 


Hot-Wire Equipment 

The hot-wire equipment used in this investigation is that designed 
by Kovasznay. (For any detailed information see ref. 20.) The equipment 
consists of functional \mits such as control \mit, calibration unit, and 
compensating amplifier, assembled together. 

The control uinit controls the heating current through the wire and 
makes it possible to measuire the resistance of the wire. 

The calibrating imit contains a direct-ciirrent potentiometer which 
makes it possible to measure the current throiigh the wire. 

The compensating amplifier is the most critical part of the hot-wire 
equipment. The amplifier in this equipment is a five-stage push-pull 
type with total gain, without compensation, of approximately 10,000. The 
compensation adds another factor of 250 making the total amplification at 

high frequencies a maximum of 2.5 X 10^. The compensation is resistance- 
capacitance type. For limiting the noise there is a five-step low-pass 
filter. 

Main units in the hot-wire equipment and the allied instruments are 
shown in the block diagram (fig. 20). 


Triggering and Timing 

From the discussion of the hot-wire behavior and of the expected 
voltage signals and their polarities ;ander various operating conditions 
in different flow regions in the shock tube, one is easily convinced 
that the hot-wire anemometer is admirably suited for timing transient 
wave phenomena. Depending upon its operating conditions, the hot-wire 
is sensitive to temperature and/or mass-flow fluctuations. The shock 
front involves a sudden jump in both temperature and mass flow. If the 
response of the hot-wire to a traveling shock front is differentiated, 
it gives a fairly sharp pulse which can be used to synchronize the tran- 
sient wave phenomena in the shock tube with the passage of the shock 
wave over the hot-wire. (See fig. 21.) If the hot-wire is operated such 
that it is predominantly mass flow sensitive, then the auxiliary arrange- 
ment shown in figure 22 is used in conjmction with the temperatiire vinit. 
This arrangement makes it possible to measure the speed of very weak 
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shocks (S = 1.01) and also strong shocks (S = 20). (For more details 
about the circuit, the scope, ajid advantages of this timing device see 
refs . 9 and 10 . ) 


Mounting of Hot-Wires 

The hot-wire anemometer has been used for two distinctly different 
purposes, namely, for actual flow stiidies and as a transducer for timing 
the wave phenomena generated in the shock tube. Thus two different types 
of probes were called for. 

The most successfiil timing hot-wire probe is shown in figure 25* The 
construction and mounting details are evident from the drawing (or see 
ref. 9)- Both tungsten (d = O.OOOI5 inch) and platinum (d = 0.0001 inch) 
wires have been successfully used. 

The probes used for act\xal flow studies are constructed as shown in 
figure 24. The mounting of these probes differs according to the con- 
figuration of the shock tube used. If the shock tube is used with the 
extreme end open, then the mounting arrangement shown in figure 25 is 
used; but, if the expansion chamber is used \inder partial vacuum or with 
the extreme end closed, then the probe arrangement shown in figure 26 is 
employed. Both platinum (d = 0.0001 inch) and tungsten (d = O.OOOI5 inch) 
wires have been used. 


Recording 

Essentially there are three different types of recording operations 
performed in this work. They are: 

(1) Shock-speed measxirements 

(2) Shadowgraph records 
(5) Oscillographic records 

If desired, they can all be used simultaneously, but for convenience 
of operation most of the time only two were recorded simviltaneously. 
Combinations recorded simultaneously were 

(a) The time the shock wave takes between two fixed stations and 

the shadowgraph 

(b) The time the shock wave takes between two fixed stations and 

the osci3JLograms 
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The arrangement and the procedure followed in each of these major cate- 
gories will be summarized below. 

Method of measuring shock speeds .- Two timing hot-wires are mounted 
at two stations some known distance apart. The output of the timing wire 
nearer to the diaphragm is fed to the triggering unit connected to the 
"start” of the chronograph. The output of the second station (which is 
downstream from the first station) is fed to the triggering unit connected 
to the "stop" of the chronograph (fig- 27). This chronograph records in 
intervals of 1 microsecond. The time taJten by the shock wave to traverse 
the distance between these two stations is recorded by the chronograph. 

Thus the speed of the corresponding shock is calculated. There are six 
available holes for mounting the timing hot-wires. The distance between 
the first hot-wire (located some I 5 feet down from the diaphragm) and 
the second one could be anywhere between 19-7 70 inches depending 

upon the configuration of the shock tube being used. 

To insure that the response of the two wires was as alike as possible, 
the wires were matched for their resistance and for the currents needed 
to heat both to the same temperature. As the timing hot-wire system has 
no bridge for measuring resitances, the main hot-wire set was used for 
matching wires. 

Shadowgraphs . - The arrangement used to take shadowgraphs of the 
traveling shock waves and the allied flow fields is reproduced in fig- 
ure 19 . The short-diiration (l microsecond) high- intensity point source 
of light is obtained by a sudden high-voltage (about 5^000 volts) dis- 
charge of a coaxially arranged set of six condensers, with a total capac- 
ity of 0.12 microfarad, across a narrow gap between two semi spherical 
electrodes (ref. 21). To trip the spark a sudden large ciirrent is passed 
through the spark coil connected to the central needle electrode. This 
generates a tiny spark which ionizes the gap and the intense discharge 
between electrodes follows. To take a shadowgraph of a traveling shock 
at a predetermined instant and location, one synchronizes the tripping of 
the spark with the help of a pickup device (timing hot-wire in this case) 
and a delay generator (fig. 19 )- The delay can be preset in intervals of 
10 microseconds. In the earlier setup, however, the required delay was 
generated by a simple re si stance -capacitance circuit which worked very 
well for low delay times (up to about 500 microseconds) but developed 
progressively higher scatter for larger delays. 

The light source is used either with a collimating lens or without 
one as the need dictates. The collimation of the light is achieved by 
an f = 4.5 and 500 -millimeter focal-length lens with the point source 
of li^t at the focal point of the lens. This gives a parallel beam of 
light approximately 4.5 inches in diameter, thus covering the full height 
of the glass section. 
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Oscillographic records.- For taking oscillograms the main hot-wire 
equipment was used. Most of the oscillograms are single-sweep records 
of the hot-wire response over a period of a few milliseconds. The prob- 
lem essentially is to bring the response of the hot-wire conveniently on 
the face of the oscilloscope. For that one has to start the horizontal 
sweep at a time somewhat ahead of the arrival of the wave front at the 
hot-wire station. This necessitates the knowledge of the approximate 
delay to be applied to the output of the triggering hot-wire. After this 
first problem of positioning the time axis, one has to calibrate the 
vertical scale such that the response signal from the hot-wire neither 
jumps off the screen nor is so squeezed that no quantitative measurement 
can be effected with any degree of accuracy. For deciding both of these 
questions the knowledge of the shock speed and the expected voltage sig- 
nal (fig. 13) helps. It will take an extremely long time to record a 
set of oscillograms if one has to calculate the jijmp in each case. After 
a while experience becomes a good guide. 

The block diagram (fig. 20 ) demonstrates the arrangement used to 
take oscillograms. For sweeping the signal across the oscilloscope once 
and only once, the output from the triggering \jnlt is fed to the oscillo- 
scope through a single- shot thyratron-tube circuit. If the ordinary 
oscilloscope (Dumont Type ^Oh- HR) is used, it is advantageous to inten- 
sify the beam temporarily by applying a flat-top piolse to the Z input 
(ref. 8 ). If the high- intensity high-voltage cathode -ray oscilloscope 
(Tektronix Type 513 D) is used, this is not necessary. Both types have 
been used during the co\irse of this work. A polaroid camera was used 
most of the time. If very fast sweep records were required (sweep time 
of approximately 100 microseconds), the high- intensity oscilloscope with 
an f = 1.5 lens and Tri X film was used. 

The operational procedure for taking oscillograms is as follows: 

The amplitude and frequency of the calibrating sine wave are adj\isted 
according to the dictates of the shock strength and operating condition 
of the hot-wire. A short preexposure is given. It Improves the qiiality 
of the Polaroid pictures and also gives a reference grid in the back- 
ground. Then the calibrating sine wave is driven across by the manual 
triggering while the camera shutter is kept open on B. After that the 
shutter is closed. Then the spot is moved toward the upper or lower 
part of the screen. This is decided purely on the basis of the polarity 
of the expected jump. If the expected jump is downward, the spot is 
moved to the upper half of the screen. The pressure across the diaphragm 
is. adjusted to the desired value. The triggering unit. Intensifying unit, 
chronograph, and so forth are all reset to receive the piilse. The shutter 
is opened and the puncturing mechanism is actuated by the person operating 
the camera. In this way he can coordinate the opening of the shutter and 
the firing of the shock tube. After the trace flashes across the screen, 
the shutter is closed. The corresponding time of traverse of the shock 
wave between two detecting stations is recorded by the chronograph. 
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Calibrating Facility 

Only the hot-wires used for actual flow studies need to be cali- 
brated. The wire is placed in the stream of air issuing out of a 
1 / 2 - inch-diameter open jet which is supplied by a vacuum-cleaner motor. 
The velocity of the jet can be varied by controlling the voltage sup- 
plied to the motor and is measured by a total-head tube mounted in about 
the middle of the jet, almost a diameter from the opening. The hot-wire 
is also mounted parallel to, and very close to, the total-head tube. The 
total head is connected to an inclined alcohol manometer. At each volt- 
age adjustment (velocity) the height of the alcohol column ir recorded. 
Also the hot-wire bridge is balanced for the same temperature by read- 
justing the current. The total pressure is given by 

-Ap = |pU^ 


or 


u = 



( 66 ) 


where Ah is the vertical height in inches of the ethyl-alcohol column. 
The constants of the wire Uq and Iq are found hy plotting (ref. 8) 


2 1 + a,u I — 

I — versus uu 

2av '' 


Typical calibration curves for platinum Wollaston (d = 0.0001 inch) and 
tungsten (d = O.OOOI 5 inch) wires are reproduced in figure 11. 

Another fimction of the calibration is to find n, the universal 
constant of the wire. This is done by the square-wave technique. The 
square-wave cvirrent is passed through the wire which distorts it. This 
is restored to a square wave by applying the needed compensation. 

* 2 

Knowing a^, M , and I at some setting of the velocity, one can cal- 
culate n from the relation 
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ANALYSIS OF EXPERIMENTAL DATA 


So far, experimental data bave been taken to examine the following 
aspects of the use of the hot-wire anemometer in shock-tube studies: 

(1) The response of a hot-wire to a sudden step-function type of 
change in flow conditions and the criterion for the proper compensation. 

(2) Ultimate time resolution of the hot-wire and the allied equipment. 

(5) Measurement of the absolute temperature jump across a shock wave 
passing over the hot-wire. 

(4) Separation of the simultaneous Jump in mass and temperature 
associated with a shock front. 

(5) Measurement of the shock speeds. 

(6) Measurement of the time of arrival at various locations inside 
the shock tube of the rarefaction front originated at the open end. 


Response to a Sudden Step-Function Change and Criterion 
for Predicting Proper Compensation 

The compensation involves usually a suitably arranged circuit which 
has differentiating properties and has increasing amplification for 
increasing frequencies, thus counteracting the integrating effects of the 
hot-wire response. Various arrangements have been used to achieve this 
end (ref. 20). 

In turbulent -flow investigations, the proper compensation is known 
from the calibration of the wire under known average conditions. However, 
when the hot-wire is operated in a shock tube, there is no such apparent 
recourse for foretelling the proper compensation. 

If a traveling shock wave passes over a hot-wire moimted in the shock 
tube, both mass flow and temperature change as a step fimction in time. 

In the previous discussion of the hot-wire response to a step function 
it has been established that the time constant may be determined by the 
operating conditions that would be assumed in a steady state sifter the 
jump. If the wire is operated at elevated temperature, after the jump 
its temperature falls; that is, its overheating ratio drops. This effect 
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is enlianced by the fact that the temperatvire of the flow behind the shock 
is higher than normal atmospheric temperature . So both of these changes 
reduce 


a^2 = 


Ry2 " ^o2 


R 


o2 


to values considerably smaller than the overheating ratio (i.e., a^.) that 

was started with. The stronger the shock, the more marked is this dif- 
ference. So the ccmpensations required are (eq. (52)) 


M2* = 



In spite of the voltage signal jump the current I is supposed to remain 
constant. Because of the finite impedance of the direct-current heating 
circuit, however. 


M2* = n - 


a^2 


(1 + 


(see section "Linearized Unsteady Response to Fluctuations"). For low 
equilibrium overheating ratios £^2 small values of c, this feedback 

correction is negligible. 

All of the quantities on the right-band side are experimentally 
(from calibration, oscillogram, etc.) determinable. So M2* can be 
computed purely from experimental data. By taking compensated oscillo- 
graphic records of the hot-wire response, one can, by trial and error, 
find the proper compensations. 

After the shock has passed and steady conditions have been estab- 
lished, one can obtain, for the flow behind the shock wave. 


l2 



2^2 
1 + a^2 



(67) 


by a similar procedure to that followed in the derivation of relation (5I) 
(refs. 4 and 10). 
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Here the wire constants Iq 2 and (^2)^ are modified from the ones 
obtained at normal atmospheric conditions: 


T 2 T 2 ke2 
I02 -lo 


h)o - § ^ «0 


( 68 ) 


j 


If the heating current I, the calibration constants n^ Iq, n.ntl Uq, 
eind the shock strength (i.e., the mass flow and temperature of the flow) 
are known, then from relation (70) the equilibrium overheating ratio 
can be calculated 6ind the theoretical value of the required compensation 
can be obtained, as 



These three different ways of getting the compensation (viz, computed, 
observed, and theoretical) have been compared with each other. Two dis- 
tinctly different approaches to data collecting were followed for this 
pvirpose : 

(a) Overheating ratios a^ (for the same wire) were kept constant 
and the shock strength S was varied. 

(b) The shock strength S was kept the same and the initial over- 
heating a^ was varied. 

For case (a) the results are reproduced in figure 28. The initial 
overheating was a^^ = 0 . 5 * The cvirrent and n remain the same. For each 
shock strength, the equilibrium overheating a^2 is different. The higher 
the shock strength, the lower is this value. Therefore for stronger shock 
strengths the compensation decreases. The computed, observed, and theo- 
retical values of the compensation agree within 10 to I5 percent. For 
high-strength shocks and constant-current operation the required compen- 
sations 6ire very small. 

The conclusions for case (b) were essentially the same as reached 
previously (ref. 8). Since the shock strength is kept constant and the 
initial overheating is varied, the current throirgh the wire is different 
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for each overheating ratio. After equilibrium conditions are established, 
Uy 2 higher for higher overheating ratios, but the corresponding cur- 

rent is higher too. The over-all effect of these two changes on 


Me 


* 


= n 


ay2 


is that the required compensation somewhat increases if one goes from an 
overheating ratio of, for example, a^ = 0.2 to ay = 0.9j keeping the 
shock strength the same for all the oscillograms. 

The unheated (or heated by very low current) wire has a time constant 
which is dependent on the mass flow (refs. 4 and 8). The relationship 
between compensation needed, for the same shock strength and the same 
wire, operated once at an elevated temperature and then with very small 
current through it, is 


M, 


Me 


ou 


1 + ay2 


( 69 ) 


where Mqu* is the time constant when the wire is operated without any 
current at velocity u (or mass flow pu). 

For weak shock strengths the temperature jirnips are very small. 
Therefore the voltage signal was low. Because of the noise level of the 
hot-wire equipment it was not possible to judge the right compensation 
from the compensated oscillograms. So the experimental points shown 
for Mq^*, the cold compensation, are for S ^ 1.25 (fig* 28). 

However, for weak shock strengths it is possible to find the right 
compensation by graphic means from the uncompensated exponential oscillo- 
grams of the hot-wire response. 

Notice that these observed values of compensations for cold opera- 
tion are somewhat lower than the corresponding values of M 2 *- The higher 
the shock strength, the greater is the mass flow. This reduces the 
required compensations. 

In conclusion, the hot-wire response to a step-function type of 
transient in flow conditions is exponential in time and can be compensated 
back to almost a step function if the compensation is set according to 
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where the equilihriiim overheating ratio a^2 is determined hy the assumed 

validity of the King’s heat-loss law (relation (29)) which is linear with 
temperature difference hut nonlinear with mass flow. 

Behind the main traveling shock wave the mass and temperature of 
the flow are constant. Tbe agreement between the expected and experi- 
mentally observed values of compensation shows that the hot-wire is 
essentially a linear instrument in its temperature response. However, 
this does not mean that the hot-wire is a linear instrument for large 
fluctmtions in general. These predicted compensations will apply for 
fluctuations in the flow behind the shock wave only if these fluctm- 
tions are smal 1. compared with the mean flow. 

It may be pointed out that if the temperature difference as it occurs 
in King's relation is large then even the temperatiare response may not be 
linear (ref. 16). 


Ultimate Time Resolution of Hot-Wire and 
Associated Electronic Equipment 

The electronic thermal-lag compensation of a hot-wire involves addi- 
tional amplification proportioneil to frequency (at higher frequencies). 

For perfect compensation this would require infinite amplification for 
infinite frequency. The total amplification in any amplifier must be 
finite; therefore, the frequency limit of the thermal-lag compensation 
is also finite. The finite frequency appears in the time domain as a 
finite rise time for the step function. Depending on particiilar applica- 
tions, the sharp rise (high-frequency response) or a noise-free operation 
may be more desirable. The increase in random-noise power is proportional 
to the area under the frequency- response curve. To limit the noise low- 
pass filters are used which limit the frequency response. (For more 
details see ref. 20.) The filter used in the present studies had a limit 
of about 70 kilocycles. Since the upper boimd of the passable frequency 
band is rather high, to resolve the rise time, one has to take oscillo- 
graphic records with very small sweep times. This puts severe demands 
on the photography of these fast single-sweep traces. As mentioned 
before, this problem was solved by using an f =1.5 lens in conjunction 
with a high- intensity oscilloscope and Tri X film. 

To start with one has to find the right compensation by trial and 
error, and then keeping the same operating conditions one has to stretch 
out the horizontal time base. Location of the jump clearly on the oscil- 
loscope face did not require much effort with the improved pickup and 
delay unit and the records could be repeated without perceptible scatter. 

A typical oscillogram is shown in figure 29- The relevant data are given 
in the figure legend. 
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The lower part of figure 29 is the response of the equipment to a 
10-kilocycle square wave passed through a duminy hot-wire with a time 
constant of 0.2 millisecond. This compensation is applied to restore 
the distorted square wave. The amplitude of the square wave is the 
same as that of the jump due to the shock wave. The oscilloscope refer- 
ence grid is calibrated for time and amplitude with a 50-kilocycle sine 
wave. The calibrating sine-wave signal was itself calibrated for fre- 
quency against the 60 cycles from the line. 

The rise time with square-wave technique is slightly smaller than 
the one observed by the response of a hot-wire to a traveling shock 
wave. The initial rise in the hot-wire trace is fairly steep, but after 
about 10 microseconds it rounds off. The same type of flattening was 
observed before (ref. 8). However, if the rise time is taken as given 
by two-thirds of the total rise, then this time is of the order of 
10 microseconds and the two techniques agree rather well. 


Absolute Temperature Measurements 

If the same hot-wire (characterized by the same Rf , Iq, (pu)q, 
and n) is operated twice with two different heating currents for shock 
waves of the same strength (i.e., the same Tq 2 P2^2)^ then under 

equilibrium conditions attained after the jump one gets 


[1 + a(T^2' - Tf)] = (Tv2* ~ To2)(a + ( 70 ) 


(I")%[l + a(T^2" - Tf)] = (t^2" " To2)(a + Bf^J 

where the primes refer to the two different overheating ratios. 


( 71 ) 


Taking ratios of these relations (for details see refs. 8 and 10 ), 
one gets 


(To2 - Ti) = ^ 


(I')^v2'l 

^2" - %) 

1 - (I")\2"( 

\ 2 ' - «l) 

»i| 

[(!' 

,■ - (r)\2" 

] I 


( 72 ) 


where 


^ _ a 

1 1 + a(T]^ - Tf ) 
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or for I' » I" and Ry 2 ' > 1 ^ 2 " 


(’’02 - ’^i) ^ 


Most of the quantities on the right-hand side of equation (72) are 
directly measurable. However, one has to determine Ry2” ^v2* 

from the oscillographic records of the hot-wire response. The equilib- 
rium resistance Ry2 is computed as follows: 

The height of the Jump due to the shock wave represents a voltage 
change ZSe across the wire. The root -me an- square value of the sine- 
wave calibrating signal is known in millivolts. The peak-to-peak value 
for the sinusoidal voltage is 2^ times the root -me an- square value. 
Therefore, the scale of the recoo^s is established. This experimental 
voltage jump in reality is smaller than what it would have been if the 
heating -current source had infinite impedance. The finite impedance of 
the heating circuit acts like a negative feedback by reducing both the 
sensitivity and the thermal lag. From the known circuit characteristics 
and the initial operating overheating ratio, one can find 

E = ^ ^ (7^) 

1 + 2a^e 


where 


+ Zg 


For cold operation 
From this, the 


0; therefore, E = 1 - e. 
corrected voltage jxmp /5>ecorrected 


can be found and 


AR = + 


^corrected 

1 


( 75 ) 


The positive sign is used when the wire gives a voltage jump corresponding 
to Increase in resistance (temperature response) and the negative sign, 
corresponding to decrease in resistance (mass flow response). Then 


Rv 2 = Rvl t ^ 


( 76 ) 
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When little heating current is used, the wire is predominantly tempera- 
ture sensitive and 


and, when the wire is operated at elevated temperatures, the wire is 
predominantly mass flow sensitive and 


Knowing ^2 > ^ 2 ^ > corresponding currents, one can find 

AT from relation ( 72 ). However, for these calculations a representative 
value had to he assumed for a, the temperature coefficient of resistiv- 
ity. The value of a was not measured for every hot-wire used. Some 
previous measurements showed that values of a given in handbooks are 
invariably higher than the values experimentally found. In these calcu- 
lations, for tungsten wires a = 0.0055 and for platinum Wollaston 
wires a = 0.005 have been used. 


When a shock wave passes over a hot-wire, the latter experiences a 
sudden temperature and mass-flow change. Its operating resistance sud- 
denly changes yielding a voltage jump. Depending upon their operating 
conditions, hot-wires are sensitive to both temperature and mass-flow 
fluctuations. If these fluctuations occur simultaneously (as in a shock 
wave), they can be quantitatively separated into temperature and mass- 
flow jumps (ref. 5)* 

By subtracting relations (70) and (7I), one gets 


^2" = ^1” + ^ 


(77) 


IV2’ = Rvi' - AR 


(78) 


Mass-Flow and Temperature Separation 



(79) 


where 


A = kgpl 



> 


(80) 


B = 
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where kg 2 is the heat conductivity of the air at equilibriiJin tempera- 
ture of the flow behind the shock wave and is assumed constant. 

Now, 


b2 2«Cpd 



(81) 


Therefore relation (79) becomes 


where 


P2^2 



. 2Io2^(Rv2' - 8v2") 




A ^ 

2(xRj^ 2(xR^ 


( 82 ) 


The calibration of the wire is usually conducted in the flow with normal 
atmospheric conditions (see section "Calibrating Facility"), and the 
constants of calibration are 


2 


2oJ^ 


(pu)o = 


2 nCpd 


(85) 


If these constants are compared with the shock-tube case (for the same 
wire) 
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or 



If the reference density^ is taken to be the atmospheric density, 

then p = and 



Also 


(84) 



(85) 


Therefore relation (82) becomes 


/P2^2 _ , 

l^r.K 1 

1 0 e -L 


V pf y 

OJ 

0 

H 

OJ 

2| 

("V2' - Rv2") 

1 

0 

4 

1 


( 86 ) 


If the shock tube is operated such that the density in the expansion 
chamber is lower them atmospheric, that is, < p^, then relation (84) 

cein be written as 



^ = U, 


^e2 




o 


MCA TN 3165 


55 


or 


P2 


Pf ^e2 


— Uo = U, 


.Pi 


Pi ^ 


and relation (86) 130000168 


^ u = ^ U ^ 

Pi ^ If Pi ° ^e2 


(!') IV2' - 


5 ( 1 ^ 2 ’ - ^2) 


I 

° ke 


( 87 ) 


Knowing the operating conditions^ the shock strength, and thus the 
expected temperature, one can find k^2/^ (ref. 15)* The rest of the 
quantities on the right-hand side of eq\iation ( 87 ) are experimentally 
determinable. Thus it is possible to calculate the mass-flow jump behind 
the primary shock wave. 

While the finite circuit correction for temperature records is com- 
paratively small, it is considerable for mass-flow oscillograms, since 
the initial overheating ratio a^ is large. This correction gets pro- 
gressively smaller if the operating temperature of the wire is so arranged 
that it is only moderately higher than the expected temperature of the 
flow behind the shock wave. In this case the overheating ratios will be 
low and jumps will be comparatively smaller too. But this technique of 
experimentation results in the difference “ ^ 2 ") ^^Tng a very 

small quantity. This small quantity, as a difference of two large quan- 
tities, is very touchy to calculate and thus the scatter in the calculated 
values of mass flows increases considerably. On the other hand, if the 
high operating overheating ratios are used, the difference " ^v2" 

(though still small ) gets to be comparatively larger than the case cited 
above. However, the evaluation of Rv2* itself is done from a large 
signal and a slightly wrong estimate of the amplitude of the actual jump 
or the calibrating signal and/or the feedback correction contributes to 
the scatter in the calculated values. The oscillographic records selected 
to calculate the mass flows, as reproduced here, used the combination of 
a temperature trace with very low current (platinum wire, d = 0.0001 inch) 
and mass flow traces with initial overheating ratio of a^ = 0.5 to 0.9« 

The computed mass-flow values are plotted against the corresponding mea- 
sured shock strength. The scatter of the measured quantities is approx- 
imately 20 percent from the theoretically expected values. 
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Besides such uncertainties as meas\iring the initial current under 
no-flow conditions (which is rarely true) and the use of the modified Iq 
and Uq, the author feels that the principal difficulty lies in the tech- 
nique of measurement and calculation. The equilibrium resistance ^ 2 ' 
due to the mass-flow response and 11 ^ 2 " temperature response 

fall closer together and they differ progressively less and less as the 
shock is stronger. This makes the personal error in setting the cali- 
brating voltage and reading the amplitude of the jump and/or of the cali- 
brating voltage rather an important factor, as all these plausible errors 

Pp 

have a strong effect on Ry 2 ' " ^v 2 " thus on — up. 

It may be added here that during the process of collecting data with 
any calibrated wire if its resistance showed a trend to increase by more 
than 0.01 ohm the wire was rejected for any quantitative mass-flow or 
temperature investigations. 

Pp 

A set of calculations of — up based on the same data as used for 

^1 

temperature calculations (fig. 50) is reproduced in figiore 51 . The 
scatter is quite evident. On the average the agreement with the theoret- 
ically expected values is fair. 


Shock-Speed Measurements 

In the section "Recording" the techniques of measuring shock speeds, 
taking shadowgraphs, and so forth have already been detailed. The chief 
interest in these measurements is the use of the hot-wire anemometer as 
a timing device which has made it possible to time shocks as weak as could 
be generated (S = 1.01) while it has been successfiilly operated at S = 20 
and it is believed that the upper limit is safely higher than this. The 
thinnest diaphragm material obtainable was of thickness 0.0007 inch. In 
the shock tube (6- by 4-inch cross section), it co\ild stand about 
5 Ib/sq in. With this diaphragm l/2-po\and-per-square-inch overpressure 
was used in the compression chamber and the time for the corresponding 
speed measurement was recorded. Because of the breaking characteristics 
of the diaphragm material at such low loadings there was some scatter 
from one shock to the other in the measurement of very low shock speeds. 

To keep the diaphragm material dry, it was stored in a box with silica 
gel. 

But for this difficulty in the mechanics of producing weak shocks, 
even weaker "shocks" could have been timed, since a hot-wire operated 
at elevated temperatures is sensitive even to normal speech. If con- 
trolled very weak "shocks" could be generated by accelerating a piston, 
the hot-wire anemometer could be a very useful instrument for investi- 
gating the phenomena. 
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During the course of this work the following interesting observation 
was made . It is rather common practice among shock- tube workers to use 
two thin sheets of diaphragm material^ if a single thicker sheet capable 
of withstanding the desired pressure is not available . When shock waves 
were being timed, a double thickness of diaphragm material was used not 
suspecting that anything unusual would happen. Quite surprisingly the 
recorded speeds were those of weak shock waves (almost the speed of 
sound). An oscillograph record, shown in figure 32, revealed that a 
very weak wave preceded the main shock about 1 millisecond at the meas- 
uring station 15 feet from the diaphragm. Apparently the double dia- 
phragm, in bursting, generated a double wave. This phenomenon consist- 
ently repeated itself, but disappeared when a single diaphragm of 
equivalent strength was used. 

The author would thus advise against the use of more than one sheet 
of diaphragm since it is doubtful if other timing devices will reveal 
the possible presence of a weak shock preceding the main front. It is 
evident that the stronger shock will eventually catch up with the weaker 
one; but, if this happens after the working section, it clearly could lead 
to erroneous results. 

Incidentally, since the first jump is undercompensated, this double 
trace is a perfect demonstration of the fact that for weaker shocks larger 
compensations are needed than for stronger shock waves if the initial 
overheating ratio is the same (fig. 28). 

For the measurement of shock speeds the thickness of the diaphragm 
was varied according to the overpressure used. For the case when the 
expansion chamber was partially evacuated, the thickness used varied 
anywhere between O.OOO3 sind O.OO5 inch. Below 100 millimeters of pres- 
sure, 0.005-inch thickness was used. It would have been advantageous if 
0.004-inch thickness were commercially available. 

The results of the measurements are shown in figure 33. These are 
average points. However, the scatter in the same shock speed measured 
more than once was surprisingly little. 


Rarefaction- Wave-Speed Measurements 

The origin and the cause of the rarefaction wave generated at the 
open end and the relations for the expected time intervals of its arrival 
at various locations inside the shock tube were discussed previously. 

For experimental check the hot-wire anemometer was used to detect the 
arrival of the first characteristic at the location of the hot-wire. 
Hot-wires were mounted at various distances from the open end and mass- 
flow and temperature oscillograms were taken. The time calibration was 
affected by the 10-kilocycle sine wave. The time between the appearance 
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of the shock wave and the rarefaction wave was measured from the oscillo- 
grams (fig. l). These oscillograms show the initial jump due to the 
primary shock wave and then some time later the arrival of the rarefac- 
tion cools the wire . This time interval between the arrival of the two 
wave fronts is the experimentally observed time. This time also can be 
theoretically calculated from relations ( 24 ) and (26): 


The agreement is good. Because of the nature of the calibration time 
base the personal error in reading it, however, the accuracy is not 
better than 3 to 5 percent. These measvirements were carried out under 
two different experimental arrangements as follows. 

(1) The hot-wire was kept fixed at one station and shocks of dif- 
ferent strengths were used and the corresponding oscillograms recorded. 
The experimental results eire reproduced in figure 34 (a). 

(2) The shock strength was kept the same and the position of the 
hot-wire was varied. The experimentally observed time of arrival of the 
rarefaction wave at these stations is shown in figure 34 (b). The 
arrangement used to mount the hot-wire for such studies is shown in 
figure 25. 

The passage of a shock wave over the hot-wire probe mounted in the 
glass section is demonstrated in figure 35- 

Hot-wires were also mounted along the ajcis but outside the open 
end (0 to 6 Inches). The time interval between the arrival of the shock 
wave and the rarefaction wave was measured (fig. 34 (b)). This can be 
done only so far as the rarefaction has not affected the central part 
of the shock wave. After this has happened, the total jirnip previously 
obtained inside the shock tube for the same shock strength falls in 
amplitude and is not constant. This is due to the alteration of the 
shock front due to the expansion wave. 

A few typical shadowgraphs taken at controlled delay settings are 
reproduced (fig. 36). The zero delay setting is at the open end of the 
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shock tube. The shock front Just emerging from the open end is quite 
straight (fig. 36(a)). The curving of the shock front is very evident 
and the central part is affected by the rarefaction about 4 inches down 
from the open end (cross section of the shock tube is 4 by 6 inches). 
The rarefaction clearly originates at the open end and its growth links 
it with the curving of the shock front (ref. 22 and fig. 36(b)). 

The vortex generated at the comers of the open end shows up very 
distinctly and proceeds rather slowly as compared with the particle 
velocity behind the central part of the shock front (fig. 36(c)). 


CONCLUDING RI34AEKS 


This experimental study has shed considerable light on the behavior 
of the hot-wire anemometer exposed to transient traveling wave fronts 
(shock waves^ contact surface and rarefaction waves) and their associated 
flow fields. As an instrument the hot-wire anemometer certainly has 
encouraging possibilities for investigation of turbulence or some such 
fast fluctuations in these transient flow fields. However, the type of 
the associated equipment of which the hot-wire will form an integral part 
may have to be decided according to the nature of the problem under study. 
There are two types of hot-wire sets in use: 

1. Constant current (heating circuit with large impedance) 

2. Constant temperature (negative feedback) (ref. 23) 

The constant -current arrangement is more commonly used in turbulence 
research than the constant -temperature system. The c ons t ant -ciir rent sys- 
tem has been used to advantage in the present investigation for study of 
the response of a hot-wire to a sudden step-function change in flow con- 
ditions, the ultimate time resolution of the equipment, measurement of 
required compensations, measurement of absolute temperatures, and for 
the detection and timing of various wave fronts generated in the shock 
tube. If a hot-wire is operated at elevated temperatures and the flow 
conditions Jump from zero to high fluctuations in mass flow (the con- 
ventional shock-tube operation), then the temperature of the wire falls. 
This reduces the sensitivity of the hot-wire to any subsequent mass 
fluctuations. However, if the constant-temperature system is used, it is 
the current which fluctuates and not the operating temperature of the 
wire. (Nevertheless, the ambient temperature of the flow will still 
Jump.) For such large negative feedback signals, the constant -temperature 
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system will be superior. Also one does not need to know the time con- 
stant as .the negative feedback suppresses it. It cannot be said, how- 
ever, that this arrangement will not have its own inherent difficulties. 


The Johns Hopkins University, 

Baltimore I 8 , Md., December 5^ 1955* 
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Figure 1.- Uncompensated and compensated response of a hot-wire to a 
traveling shock wave. Platimmi Wollaston wire, d = 0.0001 inch; 

Re = 19*12 ohms; Rw = 28.77 ohms; I = l6.24 milliamperes ; S = I.I75; 
location of hot-wire inside open end, 12.25 inches. A, compensated 
response; compensation, 0.I7 X 10“5 second. B, uncompensated response. 
C, calibrating signal; amplitude, 36 millivolts root mean square; fre- 
quency, 10 kilocycles. D, rarefaction wave. 
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Figure 2.- The (x-t) diagram of wave system produced in a shock tube. 
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Figure 3 *- Particle and wave velocities in a shock tube as a function of 
shock strength, a^^ = l,lkl feet per second; Air/Air. 
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Figure 4 .- Comparison of mass flows in region behind primary shock and 

contact -surface flow. 


MCA TN 3165 


65 





Figure 5.- Temperatures versus shock strength. Tp = 299° K; Air /Air 
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Figure 6.- Reynolds nvunber /Re versus shock strength. Expansion 
chamber at atmospheric pressure; compression chamber at higher 
pressure; d = 0.00015 inch (characteristic length); Air/Air. 
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Figure 7*- Reynolds number /Re versus shock strength. Expansion 
chamber partially evacuated; compression chamber at atmospheric 
pressure; d = 0.00015 inch (characteristic length); Air/Air. 
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Figure 8.- Variation of temperatiire throughout rarefaction wave for various 
times after diaphragm is ruptxired. Tp = 299 ° K. 


MCA TN 5165 



On 

MD 


Figure 9 »- Variation of mass flow throughout rarefaction wave for various 

times after diaphragm is ruptured. 
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Figure 10 .- Rarefaction wave generated at open end. 
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Figure 11 .- Typical calibration curve. 
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Figiire 12.- Heat loss at reduced pressures in expansion chamber. 
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(a) Flow behind primary shock wave. 

Figure 13 •- Voltage signals due to response to flows. Expansion chamber 
at atmospheric pressure. I = 0.2 centimeter; Rq = 10 ohms; tungsten 
wire, d = 0.00015 inch. 
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(b) Contact -surface flow. 
Figure 15 * - Concluded. 
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Figure l4.- "Neutral condition" in heating current. P = 1.66; 

= 15 *62 ohms. A: I = 4.58 millianrperes; jump upward; 

B: I = 5*08 milliamperes; junq) downward; 

C: I = 4.98 milliamperes; almost no jump; 

Calibrating signal: frequency, 5-kllocycle 

sine wave; an^jlitude, 0.6 millivolt root mean square. At neutral 
condition (C), opposite sensitivities for temperature and mass- 
flow junqp cancel each other. 


- 14.00 ohms. 
Rv = 14.55 ohms. 
Ry^ = l4.l8 ohms. 
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Figure 15 •- Oscillogram showing hot-wire response to shock wave and 

contact-surface flows. (Spark triggered because of response to shock 
wave.) Platinum Wollaston wire, d = 0.0001 inch; P = 75; Rp = 10.8 ohms 
a-w = 0*5* A, compensated response to shock wave. B, contact surface. 

C, calibrating signal; frequency, 1-kilocycle square wave; amplitude, 

50 millivolts.' D, firing of spark. 
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Figure I6.- Khudsen number versus shock strength. Expansion chamber at 
atmospheric pressure; compression chamber at higher pressure; 
d = 0.00015 inch (characteristic length); Air/Air. 
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Figure 17- - Khudsen number versus shock strength. Compression chamber 
at atmospheric pressure; expansion chamber partially evacuated; 
d = 0.00015 inch (characteristic length); Air/Air. 
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Figure l8.- Shock tube. 
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Figure 19- - General arrangement of equipment. 
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Figure 20.- Block diagram of hot-wire equipment 
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Figure 21 .- Timing and triggering \mit . 
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Figure 22 .- Auxiliary timing and triggering unit. 
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Figure 23 .- Timing hot-wire probe. 
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Figure 2k. - Hot-wire probe. 
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Figure 25-- Moimting of hot-wire probes. 
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Figure 26.- Probing hot-wire with vacuum seal. 
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Figure 27.- Arrangement for measuring shock speeds or setting and recording 

time delay for shadowgraphs. 
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Figure 28 .- Compensations versus shock strength. Platinum Wollaston wire^ 
d = 0.0001 inch. For n^ Iq^ and Uq see figure 11. 
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Figure 29 .- Ultimate time resolution of hot-wire and associated electronic 
equipment. Platinum Wollaston wire, d = 0.0001 inch; compensation, 

0.1T5 X 10"5 second; time constant of dummy hot-wire, 0.2 x 10"5 second 
= 19.18 ohms; Rv = 28.77 ohms; I = 16.I8 millianqperes; S = 1.175- 
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Figure 50.- Temperatures behind shock wave versus shock strength. 
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Figure 51 •- Mass flow behind shock wave versus shock strength. 
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Figure 32.- Hot-wire response to shock fronts generated by use of a 

double diaphragm^ each sheet 0.001 inch in thickness. Tungsten wire^ 
d = 0.00015 inch; = 11 . 58 ohms; R^- = 17.38 ohms; I = 27.2 milliamperes ; 

compensation, 0.24 x 10"5 second; P = 1.4. Calibrating signal; amplitude, 
^5 niilli volts root mean square; frequency, 10 kilocycles. A, weak shock; 

B, strong shock. 
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Figure 35 •- Shock wave velocities as a fvmction of shock strength (timing 
with hot-wire technique). a^_ = l,l 4 l feet per second; Air/Air. 
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(a) Location of hot-wire 12.0 inches inside from end of shock tube. 
Figiire 5^-- Detection and timing of rarefaction waves with hot-wire technique. 
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(b) Location of hot-wire varied. 
Figure 5 ^+.- Concluded. 
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L-85656 

Figure 55 •- Shock front approaching hot-wire probe. S = I. 58 . 
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’a) Shock front emerging 

(h) Location of shock 

L-85657 

(c) Location of shock 

from open end; 

from open end. 

from open end. 

S = 1.58; delay, 

1.75 inches; 

7.75 inches; 

15 X 10“° second. 

S = 2.5; delay. 

S = 2.5; delay. 

84 X 10“° second. 570 X 10“° second. 

Figure 56.- Typical shadowgraphs taken at controlled delay settings. 
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